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Abstract: Subject of this paper is ASN-Minimax (AM) double sampling plans 
by variables for a normally distributed quality characteristic with unknown stan- 
dard deviation and two-sided specification limits. Based on the estimator p* of the 
fraction defective p, which is essentially the Maximum-Likelihood (ML) estimator, 
AM-double sampling plans are calculated by using the random variables p\ and p* 
relating to the first and pooled samples, respectively. Given pi, p2, a, and /3, no 
other AM-double sampling plans based on the same estimator feature a lower max- 
imum of the average sample number (ASN) while fulfilling the classical two-point 
j- T y , condition on the corresponding operation characteristic (OC). 
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^ ! 1. Introduction 

When carrying out sampling inspection for a normally distributed character- 
s' ■ istic X ~ iV(/i, a), a > the following four cases arise: 
O. 

(i) One-sided specification limit, a known 



X 



(ii) Two-sided specification limits, a known 

(iii) One-sided specification limit, a unknown 

(iv) Two-sided specification limits, o unknown. 

In this paper, we deal with ASN-Minimax (AM) double sampling plans for 
case (iv). Let L be a lower and U an upper specification limit to X. The 
fraction defective function p(fi, a) is denned as: 

p(/i, a) := P(X <L) + P{X>U)=<$> (^A + $ (^—^] , (1) 



a J \ a 



where <£> denotes the standard normal distribution function. Note, p(fj,, a) is a 

three-dimensional function. For different levels of p, corresponding iso-p-lines 

• 11 L + U 

arise symmetrically to fi = — - — on the /z-a-plane. A figure containing 

different iso-p-lines can be found in Bruhn-Suhr and Krumbholz (1990). 

Given a large-sized lot, a single sample Xi, ...,X n , (n > 3) with 

x = ±j2 x » 5 a = - L r E(^-^) 2 , 

n ^-^ n — 1 ^-^ 

an acceptable quality level p±, a rejectable quality level P2 and levels a and j3 of 
Type-I and Type-II error, respectively, Bruhn-Suhr and Krumbholz (1990) 
develop single sampling plans based on the essentially Maximum Likelihood 
(ML) estimator 

p ' = p (X,S) = *(^-) + *(^-). (2) 



The lot is accepted within the single plan (n, k), if p* < k. 
With the help of the operation characteristic (OC) of single sampling plans, 
Vangjeli (2011) develops AM-double sampling plans A^ based on the indepen- 
dent random variables p\ and p\, which relate to the first and second samples, 
respectively. Given pi, p2, a, and f3, the AM-double sampling plan fulfills the 
classical two-points-condition on the OC and features the lowest maximum 
of the average sample number (ASN). A* is computed in a similar fashion to 
the corresponding single sampling plan (n, k) by using its one-sided approxi- 
mation AM-double sampling plan Ai, which is based on information obtained 
only from the second sample in the second stage. A double sampling plan 
consisting of two independent consecutive samples needs a larger sample size 
to fulfill the classical two-points-condition on its OC than the corresponding 
double sampling plan defined by taking into account information from both 
samples in the second stage. 

In this paper, we introduce the AM-double sampling plan A2 based on the 
random variables p* and p* Using the random variable p* which contains in- 
formation from both samples in the second stage, the OC of an arbitrary double 
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sampling plan A 2 becomes more complex than the OC of the corresponding 
double sampling plan Ai. The probability for accepting the lot after the inspec- 
tion of the first sample is analogously to Ai a single-sampling-plan-OC. Thus, 
in the next section some preliminaries regarding the single-sampling-plan-OC, 
as well as notation and definitions concerning the double sampling plan A2 are 
introduced. The increased complexity of A 2 -OC compared to Ai-OC is found in 
the probability for accepting the lot after the inspection of the second sample. 
The derivation of this probability is described in Section 3. The AM-double 
sampling plan A 2 is computed analogously to A* by using the corresponding 
one-sided approximation AM-double sampling plan A 2 . A comparison between 
A^ and A 2 is presented in Section 4. 



2. Preliminaries 

Before introducing the notation and definitions for deriving the double-sampling- 
plan-OC, we first note a well-known issue from single sampling. Let 

L( n , = P(p* < k) (3) 

denote the OC for the single plan (n, k) and let g r be the density function of 
the x 2 distribution with r degrees of freedom. 

Theorem 1: It holds that: 



with 




For the proof of Theorem 1, Bruhn-Suhr and Krumbholz (1990) use the 
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fact that for a given p (0 < p < 1) and a > 0, 



M(«r,p) := {/iGR | p{cr,n) < p} 



(5) 



is equivalent to 



M(a,p) 







otherwise, 



(6) 



with 



L-U 



(See Figure 1) 



(7) 



4 nO 




L 2 A(&>p) 



6 ^(o,p) 8 U 10 



Figure 1: Iso-p-line for p = 0.1 with /i = 5, cr = 2.431827 and M(a,p) for 
a = 1.560192 



Now, we turn our attention to the double sampling plan A2. Let X 1: ...,X ni 
be the first and X ni+1 , ...,X ni+n2 the second sample on X. Then, define the 
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following notation: 



i=l \i=l 



x, - y x,. (io) 



-2 

i=m+i 



-I ni+n.2 "T7 . "T7 

i= £ x, = 2i£i±n£?, (a) 

«i + n 2 f—' ni + n 2 

i=i 



ni+U2 

ni+1 — < 12 > 

i=i 

Definition 1: The double sampling plan by variables A 2 = 



ni ki k 2 

with ni, n 2 G N; ni,n 2 > 2; fc 2 , A; 3 G 1R + ; fci < fc 2 , is defined as follows: 

(i) Observe a first sample of size n\ and compute p\ = p(X ± , Si). 

If pi < ki, accept the lot. 
If p\ > k 2 , reject the lot. 
If h < pl< k 2 , go to (ii). 

(ii) Observe a second sample of size n 2 and compute p* = p(X, S). 

If Pp < &3, accept the lot. 
If p* > k s , reject the lot. 

The A 2 -OC is given by 

L A2 (^,a) = P M (A ± ) + P M (A 2 ) (13) 

with 

Ai = {pi < h}, A 2 = {p* p < k 3 , h<pl< k 2 }. (14) 



From (3), (4) and (14) it follows that 

P^ ff )(A 1 ) = L initkl) (fi,a). (15) 

Since P(A 2 ) '■— P^ a )(A 2 ) is more complex, we describe how to determine it 
in the next section. The A 2 -ASN is given by 

N X2 (ii,a) = n 1 + n 2 P( IJl: a)(ki <p\< k 2 ) (16) 

with 

P(n,a)(h <pl< k 2 ) = L iniM) (fi,a) - L {nukl) (n,<j). 

Remark 1: The following analogies between Ai and A 2 hold: 

(i) L\ 2 (fi,a) and N\ 2 (n,a) are not unique functions in p, but bands. 

(ii) Let the symbol * indicate the AM-double sampling plan. Denoting cf)\ 
as the one-sided AM-approximation for A^, Vangjeli (2011) shows that 
there are nonessential differences between N max (\\) and N max (cf)\) 1 . 

3. The P(A 2 ) 

Let 

P(A U 2 ) := P M (A U 2 ) = P M ( P ; < k 3 , p{ < k 2 ) (17) 

and 

P(A l 2 ) := P M (A l 2 ) = P M (p* p < k 3 , pi < h). (18) 

The probability 

P(A 2 ) := P M {A 2 ) = P M (p* p < h, h<pl< k 2 ) (19) 
can be written as 

P(A 2 ) = P(A U 2 )-P(A 2 ). (20) 

1 The examples given in the fourth section confirm this fact for N max (\2) and N max ((j)2) 
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For i = 1,2, let 



Yi-.= V^i^-^ ~N(0,1) (21) 



and 

W % := n -^S\ (22) 



a 2 



Krumbholz and Rohr (2006) have shown that the following holds: 

a y/im + n 2 ) {W l + W 2 ) + (^T 2 Y 1 -^T 1 F 2 ) 2 

O — : . 

^/(ni + n 2 - l)(ni + ra 2 ) 
Along with (21), it can be shown that 



(23) 



= o-(y x + ^n7-) + v /n 2 "o-(y 2 + v^2 ~) 

X= Z (24) 

rii + n 2 

Due to total probability decomposition and the independence of X and S 2 , 
P(A%) can be written as: 

oo oo oo oo 

P(A%) = /(/ (/ (/ PiA^W, = w u Y 1 = y u Y 2 = y 2 , W 2 = w 2 ) x 

— oo — oo 

x g n2 -i(w 2 ) dw 2 ^j $' (y 2 ) dy 2 ^j $' (y ± ) dy^j gr ni _i(wi) diui. (25) 

It holds that: 

= wi, n = yi , Y 2 = y 2 , W 2 = w 2 ) = 

= P(p(X, S) < k 3 , p(X u ft) < k 2 ). (26) 

From (6) and (24), for S < a (k 3 ), we get: 

p(X,S)<k 3 & ii(S,k 3 )<X< ti(S,k 3 ) & C 1 <Y 2 <C 2 , (27) 

where 

= (wi + W2) Mg, A: 3 ) - (gy/^T *1 + (wi + n 2 ) fx) 

1 <J v^2 



and 

(m + n 2 ) fx(S, k 3 ) - (a^nT Y x + (m + n 2 ) //) 

o 2 -j= • i/yj 

a ^/n 2 

From (23) and S < 00(^3), it follows that 

W 2 < D, (30) 

where 

(/ 7 \ \ 2 

u ^ 11,. 

rii + n 2 

Similarly, from p(Xi, Si) < k 2 , we get: 

E 1 <Y 1 < E 2 (31) 

with 



and 



Wi<F= (^) 2 (ni-l). (32) 



Setting Wi = Wi, Yi = yi, Y 2 = y 2 , W 2 = w 2 , S = S(wi, y u y 2 , w 2 ), d = 
C 1 (w 1 , yi, y 2 , w 2 ), C 2 = C 2 (w 1 , y u y 2 , w 2 ), D = D(w u y u y 2 ), E 1 = E^wi) 
and E 2 = E 2 (wi), P(A%) can be written as: 

F E 2 (w 1 ) oo 0(^1,2/1,3/2) 



P(A U 2 ) = J[ J [J ( J H(wi, yi ,y 2 ,w 2 ) x 



Ei(i«i) -00 

x g n2 -i{w 2 ) dw 2 ] $' (2/2) dy 2 ) (yi) dyi ) fc-iK) du>i, (33) 



with 
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H(w 1 ,y 1 ,y 2 ,w 2 ) = ^(C 2 (w 1 ,y 1 ,y 2 ,w 2 )) - $(C 1 (w 1 ,y 1 ,y 2 ,w 2 )). 
P(A l 2 ) is obtained by substituting k\ for k 2 in P(A 2 ). Thus, we can state: 



Theorem 2: It holds that: 

L x ^a) = L (nu kl) (^a) + P(A U 2 ) - P(A l 2 ). (34) 

4. The computation of the AM-double sampling plans 

For a given p 1 , p 2 , a and j3, the plan A^j is computed in a similar way as A*. 

rii k\ 
n 2 k 3 



ri\ k\ k 2 

We use the one-sided approximation A2 = ~ with 



= $ 1 -?lv = $ f-^=y j&s = $ ^ /:; 



where (fr 2 = I Ul ^ ^ 2 ] denotes the AM-double sampling plan in case of an 
n 2 h 

upper tolerance limit U (cf. Krumbholz and Rohr (2009)). (f) 2 is determined 
by 

(i) L^pi) > 1 - a 

(ii) ^ 2 (p 2 )</3 (35) 
(hi) N max (4>* 2 ) = min A^" max (0 2 ), 

where Z is the set of all double sampling plans <p 2 fulfilling (35) (i) and (ii). 
The AM-double sampling plan A2 is given 

(i) min L A .(<r;pi) > 1 - a 

0<CT<(7o(p) 

(ii) max L x 4a;p 2 ) < (3 (36) 

0<o-<cr o (p) 

(hi) N max {4>* 2 ) = min N max ((f) 2 ). 
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Example 1 

For L — l, U = 9, pi = 0.01, p 2 = 0.06, a = p = 0.1, we get: 

(i) (n, k) = (36, 0.02645943143) and a* = 0.082, (3* = 0.1, 

. 26 0.017577 0.035291 , 
ii A? = and N max (\;) = 32.75439. 

1 20 0.029275 



a** 


13** 


A 2 




min Lr 

a 2 


maxlr (a;p 2 ) 


0.082 


0.1 


23 
17 


0.013909 
0.026289 


0.038143 


30.45689 


0.8930783818 


0.0970618822 


0.080 


0.1 


23 
17 


0.013597 
0.026400 


0.038833 


30.72159 


0.8955304122 


0.0969993958 


0.078 


0.1 


23 
18 


0.013993 
0.026496 


0.038763 


30.99066 


0.8978607677 


0.0971362848 


0.077 


0.1 


23 
18 


0.013838 
0.026558 


0.039100 


31.12727 


0.8990880486 


0.0971046378 


0.076 


0.1 


23 
18 


0.013681 
0.026617 


0.039455 


31.26779 


0.9003201617 


0.0970742118 



where An 



23 0.013681 0.039455 
18 0.026617 



with N max (X*) = 31.26778533. 



Example 2 



For L — l, U = 9, pi = 0.01, p 2 = 0.03, a = (3 = 0.1, we get: 

(i) ( n , k) = (115, 0.0178762881) and a* = 0.085, (3* = 0.1, 

( 81 0.014029 0.021742 \ 
ii) A? = and N max {X[) = 103.5432. 

1 ; 1 I 66 0.018537 / 1 V 
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a** 


13** 


A 2 


N m ax{<t>2) 


a 2 


maxLr (a;p 2 ) 

<J 2 






72 


0.012337 


0.023495 








0.085 


0.1 








98.51959 


0.8979623972 


0.0991345737 






58 


0.017830 














72 


0.012364 


0.023535 


98.97047 






0.084 


0.1 








0.8990715849 


0.0991520345 






59 


0.017851 














72 


0.012385 


0.023569 








0.083 


0.1 


99.43030 


0.9001786758 


0.0991672779 






60 


0.017875 











where An 



72 0.012385 0.023569 
60 0.017875 



with N max (X* 2 ) = 99.43020285. 



Remark 2: Numerical investigations indicate: 

(i) The AM-double sampling plan A^j is more powerful than the AM-double 
sampling plan A* as it appears that 

N max (Xl) < N max (Xl). 

(ii) Let Ai denote the AM-double sampling plan based on the MVU esti- 
mators pi and p 2 of p(fi,a). pi and p 2 are superior over p\ and p%, 
respectively, so that 

N m ax(Xi) < N max (Xl). 

For some constellations, it could further be shown that 

N max (X* 2 ) < iNU^Ai) < N max (X{) (See Figure 2). 

(iii) The lowest N max among the AM-double sampling plans for a normally 
distributed quality characteristic with two-sided specification limits and 
unknown a would feature the plan A 2 based on the MVU estimators p\ 
and p p of p(fi,a), provided that a formula for determining the A2-OC 
would be found. 
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120 



N(p) 




Figure 2: ASN bands for A^, Ai and A£ defined by p 1 = 0.01, p 2 = 0.03 and 
a = (5 = 0.1 
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